Extremal collision sequences of particles on a line: 
optimal transmission of kinetic energy 
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The transmission of kinetic energy through chains of inelastically colliding spheres is investigated 
for the case of constant coefficient of restitution e = const and impact-velocity dependent coefficient 
e{v) for viscoelastic particles. We derive a theory for the optimal distribution of particle masses 
which maximize the energy transfer along the chain and check it numerically. We found that for 
e = const the mass distribution is a monotonous function which does not depend on the value of 
e. In contrast, for e{v) the mass distribution reveals a pronounced maximum, depending on the 
particle properties and on the chain length. The system investigated demonstrates that even for 
small and simple systems the velocity dependence of the coefficient of restitution may lead to new 
effects with respect to the same systems under the simplifying approximation e = const. 

PACS numbers: 45.50.T, 61.85 
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I. INTRODUCTION 

Chains of nonlinear interacting particles have been of 
large interest since a long time and a variety of interest- 
ing effects occurring in those systems has been described, 
such as solitons, (e.g. ||l|), energy locahzation (e.g. [g|), 
etc. In the context of granular materials chains of inelas- 
tically colliding particles have been investigated as model 
systems for shaken granular material (e.g. |g|,y), granular 
compaction g], the "inelastic collapse" (e.g. ||§,0|). The 
kinetic theory of one-dimensional granular systems has 
been addressed in j|] . 

In this paper we consider a linear chain of inelastically 
colliding particles of masses rrii, radii Ri (i — . . .n) with 
initial velocities uo = ^^ > and Vi = {i — I . . . n) at ini- 
tial positions Xi > Xj for i > j with Xi^i —Xi> Ri+i +Ri 
(Fig. ||). The masses of the first and last particles mg and 
rUn are given and we address the questions: How have the 
masses in between to be chosen to maximize the energy 
transfer from the first particle of the chain to the last 
one? If n is variable, how should n be chosen to maxi- 
mize the after-coUisional velocity w^ of the last particle. 




Sketch 



One can easily study the chains of ideally elastic 
spheres and of spheres interacting via a constant coef- 
ficient of restitution. It is much more complicated to 
deal with chains of of viscoelastic particles, which have 
an impact velocity dependent coefficient and which, as 
we show below, exhibit quite unexpected behavior. It 



has been demonstrated recently that the kinetic proper- 
ties of "thermodynamically-large" systems of viscoelastic 
particles differ significantly from those of particles inter- 
acting with constant coefhcient of restitution M. The 
system considered in this paper may serve as an example 
of a small system which properties change qualitatively 
when the viscoelastic properties of the particles are taken 
into account explicitly. 

In the present study the problem of the most efficient 
energy transmission in a chain of particles of variable 
mass is addressed. We analyze the optimal distribution 
for the particle masses and calculate the optimal size of 
the system. 



II. ELASTIC PARTICLES 

The textbook problem of elastic collisions may serve 
us to introduce the notation. Assume particle collides 
with the resting particle 1. Then after the collision the 
velocity of particle 1 is 
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(the primed variables refer to after-coUisional velocities) 
and for a chain of n-|- 1 particles of masses mo, mi . . . m„ 
one has analogously |lO| 
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For this system one finds easily that the choice mi — 
y/rrii-imi+i, {i — 2 . . .n^ 1) maximizes w'^. If we fix ttiq 
and ?Ti„, obviously the mass distribution 
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The function R^ 
has the Hmit 
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v'^/vq always increases with n and 
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i.e., if the masses of the particles are chosen according to 
(0) the kinetic energy of the first particle is completely 
transferred to the last one by a chain of infinite length. 

For the case of dissipative collisions an infinite chain 
cannot be optimal since in each collision energy is dis- 
sipated. Hence, we expect an optimum for the chain 
length for which the velocity of the last particle reaches 
its maximum. 



III. PARTICLES WITH A CONSTANT 
RESTITUTION COEFFICIENT 

According to our model the particles collide pairwise. 
This allows to use the restitution coefficient, which re- 
lates the relative velocity of colliding particles i and i + 1 
after collision to that before the collision: 
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Equation (hi) turns then into 
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where we again assume that the particle with velocity vq 
and mass mo hits a particle of mass mi at rest, which 
starts moving with the velocity v[ . Straightforward gen- 
eralization of the previous analysis for the case of the 
dissipative collisions with a constant coefficient of resti- 
tution e shows that the optimal mass distribution is iden- 
tical to that for the elastic case (pi) . This means that the 
optimal mass distribution does not depend on the dissi- 
pation if e = const. The velocity of the last particle in 
the chain reads for this case: 
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Figure shows the optimal mass distribution for different 
chain lengths n. The mass of the first particle is toq — 1 
and of the last particle ?7i„ = 0.1. 

In the next section we will consider particles which in- 
teract via a velocity dependent coefficient of restitution. 



Since the velocity of the particles varies for the particles 
of the chain we characterize the dissipation of the collid- 
ing spheres not by the coefficient of restitution itself but 
rather we define a dissipative constant b. For the case of 
a constant coefficient e it is defined as 6 = (1 — e). 
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FIG. 2. Optimal mass distribution rrii, i = 1 . . .n, for 
the case of a constant restitution coefficient e. Each of 
the lines shows the mass rm over the index i for a spec- 
ified chain length n. The masses of the first and last 
particles are fixed at mo = 1 and m„ = 0.1. 
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FIG. 3. Velocity distribution of particles in chains with 
the optimal mass distribution (given in Fig. bl) according 
to (pi). Each of the lines shows the velocity Vi over the 
index i for a specified chain length n. The dissipative 
constant is 6= (1 — e) = 5- 10~*. The last particle 
reaches its maximal velocity for chain length n* = 44 
(bold drawn). The velocity of the first particle of the 
chain is hq = 1. 

In contrast to the mass distribution the correspond- 
ing velocity distributions do depend on the value of the 



restitution coefRcient e. Figures ^ and ^ show the veloc- 
ity distribution for two different values of the dissipative 
constant, 6 = 5- 10"'* and b = 0.032. 
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FIG. 4. Same as Fig. | but for b = 0.032. The optimal 
chain length is n* = 12 
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FIG. 5. The optimal chain length n*, which gives the 
maximal transmission of energy along the chain with the 
fixed first and last masses, as a function of the dissipative 
parameter b — {1 ~ e). The line shows the prediction of 
Eq. (H), with a;o found numerically. Points refer to the 
results of a direct numerical optimization of the masses 
in the chain. 



Correspondingly, the extremal value of the R^ reads 

" 1 + e ""* 






i + XQ 



(11) 



In Fig. 1^ the dependence of the extremal n* on the resti- 
tution coefficient is shown. 



IV. VISCOELASTIC PARTICLES 

A. Collisional law for the viscoelastic particles 

It has been shown that for colliding viscoelastic spheres 
the restitution coefficient depends on the masses of the 
colliding particles and also on their relative velocity Vij 
pi| . An explicit expression for the coefficient of restitu- 
tion is given by the series |13,H 
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where Y is the Young modulus and v is the Poisson ra- 
tio. The effective mass and effective radius is defined 
a.s R'^ff = RiRj/{Ri + Rj), m'^ff — ■mimj/{mi + rrij) 
where _Rj /j and irii/j are radii and masses of the collid- 
ing particles. The constant A describing the dissipative 
properties of the spheres depends on material parameters 
(for details see |ll[). The constants Ci = 1.15344 and 
C2 = 0.79826 were obtained analytically in Ref. |13] and 
then confirmed by numerical simulations. 

For the following calculation we neglect terms O [v"^^^) 
and of higher-orders. Moreover we also assume for sim- 
plicity that all particles are of the same radius R, but 
have different masses f(l4|. We abbreviate 
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For the case of dissipative collisions the ratio Ry = 
v'^/vq does not monotonously increase with n, but rather 
it has an extremum which shifts to smaller chain lengths 
with increasing dissipative parameter h. The optimal 
value of n, which maximizes R^ reads: 



with 



^ _ log {mn/mo) 
log {xo) 
where Xq is the solution of the equation 

(l + xo) = (l-H6)a;S''/^^+^°^ 
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Thus, the collision with e = const and given dissipative 
constant b, as introduced above, corresponds (i.e. has 
equal value of e) to the viscoelastic collision with the 
same b, with unit effective mass ra'^0 — \ and unit rela- 
tive velocity Vij . 



Hence, for viscoelastic particles the velocity of the fc+1- 
rst particle after colliding with the k-th reads 
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The masses toj, , fc = 1 . . . n — 1 which maximize u^ can 
be determined numerically and the results are shown in 
Figs, y and for two different values of the dissipative 
constant b. 




40 60 

position i 



100 



FIG. 6. Optimal mass distribution m^, i = 1 . . . n, for 
the case of viscoelastic particles with the restitution co- 
efficient given by (114) with 6 = 5- 10~*. Each of the lines 
shows the mass rrii over the index i for a specified chain 
length n. The masses of the first and last particles are 
mo = 1 and m„ — 0.1. 



chain as shown in Fig. |2[ Again we find a monotonously 
decaying function for the masses. For larger chain length 
n or larger dissipation b, however, the mass distribution 
is a non-monotonous function. The according velocities 
of the particles in chains of spheres of optimal masses are 
drawn in Figs. and U. Note that the mass distribution 
and velocity distribution are related by Eq. dlq). 
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FIG. 8. Velocity distribution for viscoelastic particles 
in chains with the optimal mass distribution given in 
Fig. o. Each of the lines shows the velocity Vi over the 
index i for a specified chain length n. The dissipative con- 
stant is & = STO"*. The last particle reaches its maximal 
velocity for the chain length n* = 36 (bold drawn). The 
velocity of the first particle of the chain is uo = 1. 
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FIG. 7. The same plot as Fig. | but for 6 = 2 • 10"^. 

For small chain length or small 6, respectively, the opti- 
mal mass distribution are very close to that for the elastic 
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FIG. 9. The same plot as Fig. H but for the mass dis- 
tribution according to Fig. R (6 = 2 • 10^'^). The optimal 
chain length is n* — 20. 



B. Variational approach to the optimal 
mass-distribution 



In the following we describe an approximative theory 
of the optimal collision chain of viscoelastic particles. To 
this end we first evaluate the loss of kinetic energy in the 
chain which we divide into two parts and term as "iner- 
tial" and "viscous" losses. In our approach we treat the 
part of energy which is not transformed from the first 
particle of the chain to the last one as a "lost" energy. 
In this sense the energy is "lost" according to two mech- 
anisms: First, due to mismatch of subsequent masses, 
which causes incomplete transfer of momentum even for 
elastic collision when the masses differ (this part of the 
energy loss is called "inertial" ) . The second refers to the 
dissipative nature of collisions and, therefore, this loss is 
called "viscous" below. The inertial loss in the collision, 
attributed to the energy transfer to the j-th particle is 
thus given by the energy which remains in the i — 1-rst 
particle after the collision: 
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(17) 



For long enough chains we approximate the discrete mass 
distribution by a continuous one, m{x). This, with 
the assumption of small mass gradients, gives rrii ~ 
rrii-i + ™^' ■ 1, where we assume that particles are sepa- 
rated on a line by a unit distance. Within the continuum 
picture Ai?|^' — > ^^^ • 1, and we write for the "line- 
density" of the inertial loss, discarding high-order mass 
gradients: 
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Viscous losses describe the energy losses according to 
the inelastic properties of the material, therefore, they 
are equal to the difference of the kinetic energy of a par- 
ticle after an elastic collision (with no dissipation) and 
that after a dissipative collision. 
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Now we assume that the dissipative parameter b is small, 
so that one can keep only the linear term, expanding 



'^-^vis with respect to b. Transforming then to contin- 
uous variables and discarding terms which are products 
of b and mass gradients (which are also supposed to be 
small), yields: 



dE„ 



^^3/5^11/5^ 



dx 23/5 

Thus, the total energy loss in the entire chain reads 
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where nix = ^. As it follows from Eq. (pl|), to eval- 
uate Etot one needs the velocity distribution v{x). As a 
zero-order approximation we use an "ideal chain Ansatz" . 
This refers to a velocity distribution v{x) in an idealized 
chain, where the kinetic energy completely transforms 
through the chain, i.e. where ^m{x)v'^{x) = const = 
^rriQVQ. With mo = 1, vq = 1, so that v{x) — \/ ^Jm,[x)^ 
this Ansatz yields: 
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The mass distribution which minimizes Etot satisfies the 
Euler-equation applied to the integrand in (p2|): 
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Eq. ( p3[ ) leads to an equation for the mass distribution 
of the optimal chain, written for y{x) = l/m{x): 
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Multiplying Eq. ( p4|) by 2 (y' /y^) [y is always positive) 
we recast (24) into the form: 



dx 
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which implies the first integral of this equation: 
{y'lyf-A-2^/Hy^'^ = -c, 



(25) 



(26) 



where the constant c depends on parameter 6, the chain 
length n and initial and final masses, mo and to„. The 
form of the solution depends on the sign of this con- 
stant. If the mass distribution has an extremum at 
X — X* , such that m'{x*) — and y'{x*) = 0, the 
constant c is positive. This follows from Eq. (EQ), i.e., 
c = 4 • 22/55^1/2(2;*) > 0, since y^/'^{x*) is positive. 

The solution of the first-order equation (Bq) may be 
found straightforwardly. The general solution is some- 
what lengthy, but for the case of toq = 1 (one can always 
use the appropriate mass unit), this reads (for c > 0): 



y{x) = m{x) 
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where 
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The value of the constant c may be found from the sec- 
ond boundary condition y{n) — l/rrin, which yields a 
transcendental equation for c: 
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The last equation has to be solved numerically. Instead, 
however, we solved numerically directly the initial differ- 
ential Eq. (ph. 



Note, that some scaling properties of the solution may 
be deduced just from the form of (|24|). Namely, as it 
follows from this equation, the solution should depend 
on the reduced length variable xVb. Thus, the distri- 
bution of masses for chains with different chain length 
n and different dissipative constant 6 should coincide af- 
ter rescaling the particle numbers as i ^ v6' *, provided 
masses toq and rrin, are the same for these chains. We will 
consider the scaling properties of the mass distribution 
in more detail latter. 



Figure 10 shows the optimal mass distribution for a 
chain of length n = 40 for different damping parameters 
6. The lines display the (numerical) solution of the vari- 
ation Eq. (|24| ) whereas the points show the results of a 
numerical optimization of the chain problem. For small 
dissipation 6 both results agree. 

For larger values of 6 the solution of variational equa- 
tion (pj) deviates from the results of the numerical op- 
timization. This follows from the fact that for larger 
6 the gradients of the mass distribution are not small 
and our variational approach loses its accuracy. Note, 
however, that while the absolute values of masses in the 
mass distribution deviate from that given by variational 
approach, this still predicts well the position of the max- 
imum of the distribution. Figure O shows the same 
data as Fig. M but for larger dissipation parameter 6. 



FIG. 10. Mass distribution in chains of viscoelastic 
particles of length n — 40 with optimal mass distribu- 
tion for different values of the dissipative parameter b. 
Lines: results of the variational theory, according to Eq. 
(E4J), Points: numerical optimization (from top to bot- 
tom: • : b = 0.128, ■ : & = 0.064, ♦ : 6 = 0.032, 
A : 6 = 0.016, M: b = 0.008, T : 6 = 0.004, 
>■: b — 0.002, etc.). As previously, nii is the mass of 
the j-th particle along the chain. 
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FIG. 11. Same data and symbols as in Fig. hd but 
plotted in larger scale. 

Figure |lj displays the velocity distribution for the op- 
timal chain with the mass distribution shown in Fig. n^. 
The data given in Fig. O refer to the numerical optimiza- 
tion where Eq. (|l^) , which relates velocity and mass dis- 
tribution is used. According to the maximum in the mass 
distribution, the velocity distribution reveals for larger 6 
a pronounced minimum. 

One can give a simple physical explanation of appear- 



ance of maximum in mass distribution (and correspond- 
ingly minimum in the velocity distribution): As it is seen 
from Eq. (JlJ) the restitution coefficient increases with 
decreasing impact velocity and increasing masses of col- 
liding particles; this reduces the viscous losses. Thus 
slowing down particles, by increasing their masses in the 
inner part of the chain, leads to decrease of the viscous 
losses of the energy transfer. The larger the masses in 
the middle and the smaller their velocities, the less en- 
ergy is lost due to dissipation. On the other hand, since 
masses mg and m„ are fixed, very large masses in the 
middle of the chain will cause large mass mismatch of the 
subsequent masses and thus large inertial losses [see Eq. 
(|l7|)]. The optimal mass distribution, minimizing the to- 
tal losses compromises (dictated by b) between these two 
opposite tendencies. For the case of a constant coeffi- 
cient of restitution the relative part of the kinetic energy, 
which is lost due to dissipation does not depend on the 
impact velocity. This means that only minimization of 
the inertial losses, caused by mass gradient, may play 
role in optimization of the mass distribution. Thus only 
a monotonous mass distribution with minimal mass gra- 
dients along the chain may be observed as an optimal one 
for the case of the constant restitution coefficient. 
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FIG. 12. The velocity distribution in chains of vis- 
coelastic particles of length n = 40 with the optimal 
mass distribution according to Fig. hllfor different values 
of the dissipative constant h. Lines from top to bottom: 
6 = 2.5 • 10"'', 5 • 10"'', 0.001, 0.002, 0.004, 0.008, 0.016, 
0.032, 0.064, 0.128. The velocity distribution is obtained 
from the mass distribution (given in Fig. hlh according 
to Eq. (lid). As previously, Vi is the velocity of the i-th 
particle along the chain. 



As in the case of the constant restitution coefficient, 
the velocity of the last particle v'^ of an optimal chain 
depends on n. For short chains (with toq, ?7i„ fixed) the 
mass gradient of adjacent particles is large, hence iner- 



tia losses are large as well. For very long chains viscous 
losses become large. Hence, we expect that among the 
optimal chains exists a chain with a certain length n* 
which allows for an optimal transmission of kinetic en- 
ergy from the first particle to the last one. Fig. |l^ shows 
the velocity of the last particle for chains with optimal 
mass distribution as a function of the chain length n for 
different values of the dissipative parameter b. Natu- 
rally, as for the case of constant restitution coefficient, 
the optimal chain length n* shifts to smaller values with 
increasing dissipative constant b. 
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FIG. 13. Velocity of the last particle w„ for chains of 
viscoelastic particles with optimal mass distribution over 
the chain length n for different values of b. As in Fig. [13, 
the velocity distribution was obtained from the mass dis- 
tribution according to Eq. (Ilq), and lines from top to 
bottom correspond to 6 = 2.5 • 10"", 5 ■ 10"*, 0.001, 
0.002, 0.004, 0.008, 0.016, 0.032, 0.064, 0.128. Note that 
with increasing dissipative constant b the maximum of 
v„{n), which corresponds to the optimal chain length n* 
shifts to smaller values of n, which means naturally, that 
optimal chains are shorter for larger dissipation. 

Having the mass distribution and the velocity distri- 
bution obtained from the numerical optimization one can 
check directly the validity of the " ideal chain Ansatz" , 
v{x) = 1/ y^m{x), used in the variational approach. In 
Fig. n^ we compare v{x) obtained by optimization with 
that from the Ansatz. As it is seen from the figure, the 
ideal chain Ansatz occurs to be rather accurate for small 
dissipation parameter 6, and for the initial part of the 
chain. It demonstrates, however, noticeable deviations 
from the optimization data for larger b, especially at the 
end of the chain, i.e. for i ra n. This is not surprising 
since it uses an assumption of complete transmission of 
energy, which is definitely poor for the very end of the 
chain. On the other hand, as it follows from Figs. M 
and nil this Ansatz yields rather accurate results, when 
applied to the mass distribution problem. The possible 



explanation for this follows from the boundary condi- 
tion for the mass distribution at the end of the chain, 
m{i ^ n) = rUn- This imposes the correct behavior of 
the mass distribution at this part of the chain and partly 
compensates the inaccuracy of the velocity distribution 
which develops mainly at the chain end (see Fig. 04). 
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FIG. 14. The velocity distribution in chains of vis- 
coelastic particles of length n = 40 with the optimal mass 
distribution according to Fig. hol Lines give the veloc- 
ity distribution for the ideal chain Ansatz, Vi = l/^/mi 
(with masses taken from the optimization data), points 
show the numerical optimization data for b = 0.001 (top) 
and b — 0.008 (bottom). Note that for these values of the 
dissipative parameter b the variational theory gives very 
accurate description for the optimal mass distribution 
(see Fig. 0) 



C. Scaling laws for the optimal mass-distribution 

Now we analyze how the maximal mass m* = m(x*) 
(the mass of the heaviest sphere located at a; = x* ) in the 
optimal mass distribution depends on the chain length n 
and the dissipative parameter b. We show that there 
exists a simple scaling relation between these values. 

We start from Eq. (g7|) for the optimal mass distribu- 
tion 
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with c and (p defined by Eqs. ( p8| ) and (^9|). The condi- 
tion for the optimal mass 
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the maximal mass m* and the constant c 
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c = 2^/^b/V^. (33) 

This allows to write the boundary condition for m{x) at 
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Simple analysis shows that (^ < if the optimal distribu- 
tion has a maximum (this follows from the form of the 
solution (pfl) and the requirement that mix) increases at 
a; = 0). Thus, one obtains from Eqs. (ES) and (|33 
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Using again Eq. (33) for the constant c we recast Eq. 
(pq) into the final form: 
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This scaling relation expresses the product nv6 in terms 
of the maximal mass m* . For the case of a strongly pro- 
nounced maximum in the optimal mass distribution, i.e., 
when mQJra* <C 1 and mn/m* <ti 1, one can expand 
the arccos(a:) in ( |37| ) to obtain a linear scaling relation 
between {m*)^'^ and ny/b: 



with 
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In Fig. Oa we compare the analytical relation i pT\ ) and 
its linear approximation (Bq) with the results for to*, 
following from the numerical optimization for the mass 
distribution for different chain lengths and different dissi- 
pative constants. As one can see from Fig. n^, the results 
of the analytical theory and of the numerical optimiza- 
tion agree well, except for large dissipation values. We 



would like to stress that there are no fitting parameters 
used. 




FIG. 15. Shows n\fh as a function of m* for the chain 
of viscoelastic particles with the optimal mass distribu- 
tion. Here m* is the mass of the heaviest particle in the 
chain, n is a chain length and h is the dissipative param- 
eter. In the figure we plotted nyh over m* for about 
3000 different combinations of 6 and n {n — 2 . . . 300, 
6 = 0.0001 . . . 0.256) including all data presented in Figs. 
pi M, hd and uA Without any adjustable parameters the 
data from the numerical optimization of chains agrees 
well with the analytical expressions Eq. (p7|), given by 
the dashed line. The linear approximation for the scaling 
relation, Eq. (BSl) is shown by the dotted line. 

Using the optimal mass distribution, Eq. (^), one can 
compute the total energy loss in the chain, as given by 

Eq. My. 



T^ "^^ , n r I V^l - C0S2 ipn \/^ " COS^ (f 

t^tot — -^ + ^V C < 

2 I COSyJn COSip 

-f-2v^{arcsin [cos</?„] — arcsin [cos (^] } , 



(41) 



where (/?„ = nyc/2 + ip. According to Eq. (30) one 
obtains 



mo 



m„ 



34^52 

24/5^2 



■ COS if 



■ COS ifin , 



(42) 
(43) 



which allows to express all trigonometric functions in Eq. 
( p] ) in terms of mo and m„, yielding 



Ef. 



/22/5& 

Z172 



/ 22/56 



en 



(44) 



and finally, taking into account Eq. (33) for c, we arrive 
at the relation for the total losses 



£tot(n,6) = 2^/5^6 




(45) 



/mo Vm* 2^/5 Vm* 

Using the approximation for the maximal mass 

m* « ( nvbj-p + qj-p \ , 



(46) 



which follows from Eq. (^8|), one obtains an explicit ap- 
proximate relation for the total losses and, thus, for the 
final velocity 



,,'2 



Etot(n,b) 

m„ m„ 



(47) 



in terms of of the chain length and the dissipation con- 
stant h. Unfortunately, due to the fact that chains with 
optimal lengths obviously do not have a maximum in 
their mass distribution, one cannot use the previous re- 
lations to estimate the optimal chain length for a given 
dissipation constant 6, since these relations hold true only 
for chains which do have a maximum. 

Note that since the maximal mass m* depends only 
on the product n\h, the expression in curled brackets in 
the right-hand side of Eq. (^) also depends only on this 
combination. This suggests the following scaling rela- 
tions for the final velocity for the chains with fixed n^/h: 



,,'2 



-dVb 
- d'/n , 



(48) 



where we take into account that mo = 1, wo = 1, and 
where d and d' are some constants which are defined by 
the particular value of nyb. 



V. CONCLUSION 

We investigated analytically and numerically the 
transmission of kinetic energy through one-dimensional 
chains of inelastically colliding spheres, where the first 
and the last mass is fixed. For the case of a constant 
coefficient of restitution we found that in the chain with 
optimal energy transmission the mass of each particle is 
given by the geometric average of its neighbors, i.e. the 
distribution of the masses of the spheres is a monotonous, 
exponentially decreasing function. This function is in- 
dependent on the coefficient of restitution e where the 
special case of elastically colliding particles (e = 1) is in- 
cluded. We derived an expression for the chain length n* 
which leads for a given e to the optimal energy transfer 
(provided the masses in between the first and last mass 
have been chosen properly). 



The situation changes quahtatively if we assume that 
the chain consists of viscoelastic spheres for which the 
coefficient of restitution depends on the impact velocity. 
Here, the optimal mass distribution which leads to max- 
imum energy transfer is not necessarily a monotonous 
function. Depending on the chain length n and on the 
material parameters of the spheres it may reveal a pro- 
nounced maximum. The part of the kinetic energy of 
the first particle, which has not been transfered to the 
last one, we consider as losses of energy. These losses 
have been characterized as losses according to incom- 
plete transfer of momentum due to mass mismatch of 
the particles (inertia losses) and losses due to the dissi- 
pative nature of particle coUisions (viscous losses). We 
develop a theory which describes the total energy losses 
along the chain, so that the optimal mass distribution, 
minimizing the losses, may be obtained as a solution of 
a variational equation. We find a general solution to this 
nonlinear second-order differential equation. Implication 
of the boundary conditions yields, however, a transcen- 
dental equation, which one needs to solve numerically (in 
practice, we solve numerically the initial differential equa- 
tion). We observed that our variational theory agrees 
well with the results of the numerical optimization for 
the mass distribution, provided the dissipative material 
parameter is not too large. We also performed a direct 
verification of the basic approximation used in our vari- 
ational approach. 

From the exact solution of the variational equation we 
obtained an analytical expression which relates the heav- 
iest mass in the mass distribution to the chain length 
and the dissipation constant. We found that this analyt- 
ical expression, having no fitting parameters, is in good 
agreement with the numerical data. Using the exact so- 
lution for the optimal mass distribution we also found an 
expression for the total energy losses. This allowed to 
obtain scaling relations which show how the velocity of 
the last particle in the chain scales with the length of the 
chain n and with the dissipation constant b, for the chain 
with the value of nVb fixed. 

It has been demonstrated before that for the case of 
"thermodynamically-large" granular systems the impact- 
velocity dependence of the restitution coefficient, as it 
is given for viscoelastic particles, may lead to qualita- 
tively different behavior as compared to systems with a 
constant restitution coefficient, e.g. 0,0,^3]. The system 
investigated here may serve as an example of the major 
influence of the velocity dependence of the restitution co- 
efficient even for relatively small ( "lab-scale" ) and simple 
systems. Therefore, in general, the assumption of a con- 
stant coefficient of restitution is an approximation which 
justification cannot be assumed a priori but has to be 
checked for each particular application. 
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